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Abstract 

In the present article we show that the Skyrme-Faddeev model 
possesses nonlinear wave solutions, which can be expressed in terms of 
elliptic functions. The Whitham averaging method has been exploited 
in order to describe slow deformation of periodic wave states, leading 
to a quasidincar system. The reduction to general hydrodynamic sys- 
tems have been considered and it is compared with other integrable 
reductions of the system. 
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1 Introduction 



The main aim of the present paper is to study periodic and multi-periodic 
solutions of the so-called Skyrme-Faddeev model. In the recent years a 
special interest was deserved by the 3D static nonlinear Skyrme - Faddeev 
a- model for the field : M 3 — > § 2 , the total free energy of which is given by 



Ssf = I d 3 x 



H ik = d>- [did> x d k 4>] , , . 
0-0 = 1, { - } 



where p is a real positive constant and Hi k is an antisymmetric tensor field, 
the so-called Mermin - Ho vorticity [1, 2, 3] in Matter Physics, expressing 
the non irrotational properties of the multi-component superfluid. On the 
other hand, after a long work [4, 5], it was proved that (1.1) can be seen as a 
special subcase, among many others [6], for the background (classical) field 
of the quantum pure SU (2)-Yang- Mills theory in the infrared limit, thus 
describing a self-consistent "mesonic" field in the context of the Nuclear 
Physics. 

The main interest is to look for localized finite energy solutions, for which 
one imposes a constant value at spatial infinity, for instance limix^^ = 
(0,0,1), compactifying to § 3 the space domain. Then, from the homo- 
topy group theory result 1:3 (§ 2 ) = Z [7], one can conclude that all such 
solutions (hopfions) are labelled by the so-called Hopf index Q into sepa- 
rated sectors. The Hopf index can be computed analytically, since Hi k in 
the expression (1.1) is a closed 2-form on § 3 , so derivable from a 1-form 

a k = - E -f* f„ 3 ( x *~ x '^ Hl i ( *'l d 3 x'. Thus one has [8, 9, 10] 

Q = ~7~7r~n [ d 3 xsikiaid k ai, 



16tt 2 . 

expressing the linking number of the pre-images of two independent points 
on the target § 2 , then Q = for spherical solutions, 1 for toroidal shaped 
vortices. Only approximated analytical solutions are known [11]. Many 
times tangled hopfions have been confirmed by numerical studies [12, 13, 
14, 15, 16] , which have produced a comprehensive analysis of solitons with 
1 < Q < 16, proving the existence of local energy minima with knotted 
structure. On the other hand, global analytical considerations [8, 9, 10, 17] 
have shown that the energy of such a knot is bounded from below by 
S S f > C it 2 p IQI 3 / 4 , where C is a constant evaluated numerically [16]. 
The main consequence of the above bound is that hopfions of higher topo- 
logical charge can be broken only by adding an extra energy contribution for 
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their disentanglement. However, also space extended structures were con- 
sidered in [18] and in [19]. In the latter it was shown that periodic collection 
of localized objects may condensate in order to form periodic structure in 
the space. However, as pointed out in [19] the appearance of extended multi- 
sheeted structures, similar in some extent to the stripes found in [18], may be 
energetically more favorable. Thus, the quest for periodic (possibly exact) 
solutions for the Skyrme-Faddeev model becomes more interesting. On the 
other hand, in a series of papers [20], it was shown that, by adding certain 
differential constraints to the equation of motion of the Lorentz invariant 
version of (1.1), one may obtain completely integrable sub systems, with in- 
finitely many local conservation laws. This was another hint in the direction 
of integrable reductions. On the experience of the hydrodynamic reductions 
we tried different approaches, which lead us to identify 1) quasi-periodic 
solutions described in terms of elliptic integrals, 2) the average motion of 
periodic waves accordingly to Whitham's method, 3) constraints on the 
solutions expressed in terms of several Riemann invariants. The paper is 
organized into an Introduction and three further Sections, concerning the 
previously listed topics and final considerations are included in Conclusions. 



2 Periodic Solution 

First we consider the Skyrme-Faddeev model in the 4-dimensional relativis- 
tic space time [9], defined on the space-time M = { far, . . . , x 3 ) } endowed 
with the (pseudo)-Riemannian metric diag (g^) = (+,—,—,—), given by the 
Lagrangian density 

C = 32^ ' d ^ ~ \ ^ X • ^ X ^) - « (1 - ^ » > 

(2-1) 

where A = -f > is the scaling parameter describing the breaking of the 
conformal symmetry and K is a Lagrangian multiplier implementing the 
constraint cf) G 5^. The model (1.1) is obtained by setting do<p = 0. 

It is well known that the geometric constraint <j> ■ 4> = 1 can be realized in 
several ways, but here it seems useful to introduce the polar representation 

4> = (sin O sin $, sin 6 cos <£, cos , (2.2) 
where O and $ are suitable function on the variables far,..., a; 3 ) to be 
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determined. The Lagrangian (2.1) becomes 



(2.3) 

Now, it is well known [11] that the equations of motion given by (1.1) admit 
harmonic plane wave solutions of the of the form 

/ 2 A cos (piX*) 2A sin (piX 1 ) 1 - A 2 \ 
Ct)= \ A 2 + l ' A 2 + 1 ' A? + T ) ' (2 ' 4) 

where ^4 parameterizes the amplitude of the third component and remark- 
ably the dispersion law is given by pip 1 = . Actually, by using the symme- 
try group of the model one can find a 12 parametric family (4 translations, 
4 boosts/rotations and 3 gauge transformations) of solutions. In particu- 
lar, the axis of precession can be arbitrarily fixed by gauge transformations. 
Nevertheless, they cannot be superimposed, because of the nonlinearity char- 
acter of the equations of motion. To have a visualization of (2.4) can think 
to a periodic assembly of vectors, whose wave front are orthogonal to the 
direction p = (pi,P2,P3), and maintaining constant the projection fa. The 
configuration is similar to a cholesteric liquid crystal. Moreover the energy 
density of the configuration is constant in the whole space, being equal to 
£ = 47r 2(i_^2)2 • Notice how solutions with smaller third component are 
more energetic. 

Looking for the simplest generalization of (2.4), one assumes that 

8 = G [9] , $ = $ [9] + 9, where 9 = a^, 9 = p^af (2.5) 

in which one distinguishes 9 as the phase from the pseudo-phase 9. From 
a different point of view, we are looking for invariant solutions under a 8- 
parametric family of 2-dimensional Abelian sub-algebra of the translation 
symmetry group given by 

{3 3 3 3 1 

v a = ^2 ait ° ~ a ° X) **' V P = X] P ito - Aj ti L (2.6) 
i=l i=l i=l i=l J 

where tj's are the generators of the translations. Actually, by using the 
adjoint action of the space-time rotational subgroup, we can conjugate each 
of the above subalgebras to exactly one representative sub-algebra of the 
form (2.6) belonging to a 3-parametric sub-family. However, it is more easy 
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to deal with all components for notational homogeneity. Thus the equations 
of motion reduce to the announced 3-parametric family 



2B 3 - -£sin 2 6 



2B 3 sin 2 6 <5>oe + @e sin 26 (2B 3 <f> e + B 2 ) = 0, 



sin 2G l^Be 2 e + B 3 $> 2 e + B 2 § 6 + #i J2.7) 



(2.« 



where B 1 = -p^p", B 2 = -2a^ B 3 = -a^a" and B = B\- AB X B 3 . 

Of course these equations provide the above linear solutions (2.4), setting 
B 3 = and B 2 ± 0, then 6 = 2 arctan^, $ = $ - %Q and Pi = B ^~^ . 
On the other hand, for B 3 ^ 0, the solution is given by 6 = 2arctanA and 

<!> = $0 — B% 2B^ ®-> w ^ n an analogous expression for the p^s. 

To deal with the general situation one uses the expression of the energy- 
stress tensor T^ v = (©"d^ + <^<%J C P - g^C P . Substituting in it the 
ansatz (2.2), only derivatives with respect to 9 survive. Thus, the further 
vanishing divergence is equivalent to take dg over a quantity obtained by 
contracting £^ = T^ v a y , corresponding to total conserved quantities for the 
wave, seen as a function of the phase. Their expressions are 



-1 



32^ 2 



{B 3 a e 2 e + sin 2 6 [2a • (3 /3 + (#i - 2/3 2 ) a 



+B 3 (2p + a Q ^e) <S>e 



XB 



ao©6 



32vr 2 



B 3 aiQ 2 e + sin 2 6 



B 2 Pi - B iai + B 3 (2/3; + a&g) $ e 



XB 



(2.9) 



These equations can be used to find an expression of 6g and <&q. Precisely, 
assuming B 3 ^ one finds 



R 2 



8B 3 (Si sin 2 9 + U 3 ) - 2B\ (sin 2 6 + E7f esc 2 6) 
B 3 (8B 3 - XBsin 2 6) 
B 2 (1/2 esc 2 (6) + 1) 
2B 3 ' 



(2.10) 
(2.11) 



where the U^s are two constants completely defining the quantities in (2.9) 



by the expressions £ M = U 3 a,n + 



B 2 u 2 



Thus, one has re- 



*m 1 2 

duced the problem to the quadratures, introducing only two new integration 
constants besides {U 2 ,U 3 }, which determine the amplitudes of the phases. 
Finally similar conclusions can be obtained in the case B 3 = 0. Despite of 
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its involved expression, equation (2.10) can be set in algebraic form by the 
transformation 

6 = arcsin v / ^, (2.12) 



forcing to be < ip < 1 and to be satisfied the equation 

64(V-l)(^-Ai)(V-A 2 ) 



\ 2 B4>i (ipi - V) 



(2.13) 



where one has defined the constants Ai j2 = 2B:iUaziz V 4: B 3 u 3 B U^_^ re l a ted 1- 
1 to the values of the integrals of motion, and ip\ = jjg-- Assuming Ai 
to be real and setting < At < A2 < 1 for instance, by a continuous 
variation of ipi one obtains different behaviors of oscillation amplitudes for 
V>: A 2 < ip < 1 for ipx < 0, A x < tp < A 2 for < Vi < A 2 , A x < ip < 1 
for ipi = A2, A 2 < ip < 1 for A2 < ipi- So for all choices of ipi there is 
only one oscillating solution, bounded between two of the three zeros of the 
numerator in (2.13), even if real unbounded solution may appear or complex 
ones (see Figure 1). 

Analytically, equation (2.13) can be integrated in terms of incomplete 
elliptic integrals of the third kind. Precisely, by introducing a parametric 
variable Z one obtains the parametric form 



wo 




£A 2 V>! (V>! - A,) 2 
4 V (Ai - 1) {A 2 - ipi 



n 



1) (Ai - A 2 ) 



^-A 2 ' 1 (Ai - 1) (Vi - A 2 ) 



A 2 ^i sin 2 Z + Ai (-01 cos 2 Z - A 2 
Ai sin 2 Z + A 2 cos 2 Z + tpi 



(2.14) 



which can be expressed in terms of Weierstrass V function. Furthermore, 
from (2.11) the function <3? can be expressed again in terms of incomplete 
elliptic integrals, namely 




WW) 



+ 



+ $0 



2^i (Ax-Vi) 2 (S 1 AVi + 2) n /A 2 



(A!-1)(A 2 -Vi) 



Ai 



A 2 -Vi 



(2.15) 

(A x - A 2 ) (Vi - 1) \ 



A 2 ^i (Ai - V1) 2 



-n 



(Ai - A 2 )-0i 

Ax (A x - 1) (A 2 - VAi (</>i - A 2 ) 
where s\ = signi? 2 , s 2 = sign U2 



(A x - 1) (^ - A 2 ) 



(Ax - A 2 ) (Vi - 1) 



(Ax - 1) - A 2 ) 
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In this parametric form it is evident that the phases of the spinorial 
field and the phase 6 do not have the same periodicity. Thus the solution is 
generically quasi periodic and only for very special choices of the parameters 
true periodic solutions appear (see Figure 4 ). Then it is convenient to adopt, 
as we will show in the next section, a method which allows to describe 
solutions with a minimal set of parameters, concerning only the periodicity 
in the phase. Here notice only that the length-wave can be made very large 
when A 2 — > 1 and ^1 — >• 00. 

Then, an important observation occurs now: one has found a special 
2-dimensional reduction of the Skyrme-Faddeev model, which is completely 
integrable, and one may wonder if this is not in the class described in [20]. 
If one performs the transformation of the polar representation of the field <fi 
into the stereographic projection 

w = ztan ( ^- ] exp (— i<&) ,w* — > — itan ( — ] exp , (2.16) 



2 r 

the constraint imposed by the authors in [20] is expressed by 

d^w d^w = 0. (2.17) 

One easily verifies that it is satisfied by the harmonic wave solution (2.4). 
On the other hand, if one replaces in (2.17) the reduction in (2.5), depending 
on phase and pseudo-phase, one obtains the relation 

sin 2 (6) (B 3 $ e + B 2 ) + B x ) - B 3 @j + 

i sin(0) (25 3 $ e + £ 2 ) e = (2.18) 

from which one sees that separately real and imaginary parts have to vanish. 
Replacing the condition (2.11) into the last equation, one obtains 

2 1, , , x ^(sin^O)-?/! csc 2 (0)) 
-B s @ 2 e + B l sin 2 (0) - iB 2 U 2 e e csc(G) - AO — = 0, 

4-D3 

(2-19) 

saying that, excluding constant solutions, or B 2 or U 2 have to be zero. 
But compatibility with equation (2.10) implies the constancy of in both 
cases. So, in conclusion, the solutions found above are out of the sub-sector 
described by the constraint (2.17). 

In the next Section we average the periodic solutions of the Skyrme- 
Faddeev model by the Whitham approach. 
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3 The Whitham averaging method 



The Whitham approach was developed for any multi-dimensional system 
possessing a Lagrangian formulation (see details in [2 ]). Nevertheless, only 
some year later, the first averaging on a multi-dimensional example, the well 
known three dimensional Kadomtsev-Petviashvili equation, was provided by 
E. Infeld [22]. Now, more than thirty years later, we present the second 
multi-dimensional example, i.e. we are considering the averaging of the four 
dimensional Skyrme-Faddeev model. As we have seen above, this nonlinear 
system is determined by a Lagrangian (see (2.1) and (2.3)) and possesses 
a multi-parametric family of periodic solutions (2.13) (see also (2.10) and 
(2.11)). In such a case, following a heuristic approach, one can introduce the 
averaged Lagrangian L(j, uj, (3, k) in the new dynamical variables 7, uj, (3, k, 
which correspond to the derivatives with respect to space time variables of 
phase 9 and pseudo-phase 9, now not necessarily linear as in (2.5). It means 
that ui = —6 x o,k l = 9 x i and 7 = — 9 x o,j3 % = 9 x t. Thus, one immediately 
derives the four-dimensional quasilinear system 

d x oL u = d x iL kl , d x oL~f = d^Lpi, (3.1) 

with the compatibility conditions 

*£o + LJ xi = 0, , k xJ = kl i + j, (3.2) 

^0+7^ = 0, , 

The averaged Lagrangian density £(7, uj, (3, k) can be obtained in two steps. 
Formally, one replaces the family of periodic solutions (2.10) and (2.11) into 
the Lagrangian 2.3, obtaining 

C p = sin 2 (6) (-^X - BxB^j Q 2 e + B 3 *i + Me + B^j + B 3 Q 2 e , 

which is a function only on 9 thus, performing an integration over a fi- 
nite space-like region, contributions from ^-independent coordinates are just 
time-independent finite multiplicative factors. Then, on a period of the 
wave, one leads to the Lagrangian 

L= j> C p d9, 

where, generalizing the standard Whitham approach, we need to introduce 
two natural normalizations (or constraints) 

j> d9 = 2tt, < $ >= j> $d9 = 2vrm, (3.3) 



S 



where the integer "to" is the number of rotations of the vector </> around a 
value determined by a given pseudo-phase 9. This situation is very similar to 
to the spin wave configurations called cyclon and extra-cyclon in multiferroic 
materials [23]. Then the corresponding averaged Lagrangian density is given 
by 




(3.4) 



where we introduced the function 




w=iz*\r -^w-^w-^. (3.5) 



Thus, one immediately can check that two equations La 1 = and La 2 
(see [21]) coincide with normalizations (3.3), while the Euler-Lagrange 
equations lead to four dimensional quasilinear system of the first order (3.2). 

Remark: The two normalization conditions (3.3) formally allow to ex- 
clude A\ and A 2 from the above construction. This means that solving 
(3.3), one can derive Ax (7, uj, ft, k) and ^.2(7, to, f3, k). Thus, quasilinear sys- 
tem (3.2) contains first order derivatives with respect to x l of eight unknown 
functions (7, w,/?,k) only. However, one can use the five roots 0,1, A\,A 2 
and ipi for parametrization of periodic solution (2.13). Nevertheless, the 
function TV (7, uj, ft, k) as well as the averaged Lagrangian cannot be ex- 
pressed via complete elliptic integrals of the first and second kind only. For 
this reason, all partial derivatives L w , L k i, L 7 , Lpi contain also elliptic in- 
tegrals of the third kind, thus all expressions in quasilinear system (3.2) 
became too complicated to be presented explicitly in this paper. 



4 The Method of Hydrodynamic Reductions 

In comparison with the approach considered in the previous Section, we 
describe a special class of solutions of the Skyrme-Faddeev system, which 
can be obtained by the method of hydrodynamic reductions (see [24]). The 
Lagrangian density C p depends on 0,0^ and <$> u . Actually, the standard 
method of the hydrodynamic reductions is applicable to a Lagrangian den- 
sity which depends on derivatives M and § u only (see [25]). Nevertheless, 
even in this more complicated case, the method of hydrodynamic reductions 
can be utilized. Indeed, the Euler-Lagrange equations 

d M 0> = - sin(20)<l\<l> + ^sin© • <l^[sin 0(0^<T - 9"*")], 
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0^sin(20) + sin 2 G[d^ + -Q v d^@ v - = 

contain the field variable 0, its first and second derivatives 0^,0^, while 
the field variable <& is involved only by its first and second derivatives 
^n^iiv- According to the method of hydrodynamic reductions we intro- 
duce N Riemann invariants r l (x,t,y, z), which satisfy simultaneously the 
three commuting diagonal hydrodynamic type systems 

4 = /ArK, r { y = rf{r)rl r\ = C*(r)r{, 

where new auxiliary field variables u v = & u and ui M = M depend on these 
Riemann invariants only. For successful application of the hydrodynamic 
reduction method, we should select all parts of these Euler-Lagrange equa- 
tions containing just first and second derivatives ^^^^ and annihilating 
the coefficients of sin©. This yields the constraints 

= o, (\ + ^e„@A = o (4.i) 

implying the nonlinear system 

= o, d^+^e u d^e u -$ u &) = o, ^ v d fl (Q^ v -e v ^) = o. 

The latter is equivalent to the quasilinear system of the first order 

d^vf = 0, - ^w^u^ - vPw") = 0, Ur] d^vP - wV) = 0, 

(4.2) 

so that (4.1) reads 

upvf = 0, (l - ^WfivA = 0. (4.3) 

The last relationship contains two admissible constraints 

u^u" = 0, 1 - -wpw? = 0. (4.4) 

Thus, the method of hydrodynamic reductions is applicable for (4.2) if it is 
equipped by the constraints 

u^vf = 0, = (4.5) 
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or by the two other constraints 

2 

Utl w^ = 0, w a w a = -. (4.6) 
A 

Actually, only the second choice exists for the Skyrme-Faddeev model, while 
the first one also is possible for A = 0, that is for the usual 0(3) non- 
linear ex-model. Furthermore, let us notice the similarity of the previous 
constraints on the existence of hydrodynamical reductions for the Skyrme- 
Faddeev model with the constraint (2.17). 

5 Conclusions 

We have found exact analytic quasi-periodic spin waves for the Skyrme- 
Faddeev model. We have studied in detail their main features. In particular 
they are determined in terms of elliptic integrals of third kind. Assum- 
ing that there exist solutions which are periodic, we have shown that the 
Lagrangian can be averaged by the Whitham method. This provides a La- 
grangian for a set of parameters, describing the evolution of periodic waves 
in terms of a quasilinear system in partial derivatives of the first order. Fi- 
nally, noticing that in a general such a system like (3.2) is non-integrable. 
But, one can use the so called "method of hydrodynamic reductions", which 
allows to extract particular reductions, whose solutions can be parameter- 
ized by a set of arbitrary functions. This can be achieved only if some special 
constraints (4.5) or (4.6) are satisfied. However, the analysis of those con- 
straints deserve several technical complexities and it is not done here. Only, 
we remark that the constraints are a further restriction on the subsector of 
the solution space described in [20]. 
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Figure 1: The graphic for the inverse square root of (2.13) for the family 
of parameters B = 1, A± = .1,A% = .8 and —.45 < tfti < 1.55 with steps of 
0.1. Colors run accordingly from red to violet. Only one bounded periodic 
solutions exists for any set of parameters. The degenerate case ip\ = A2 = .8 
is not considered, but it corresponds to the confluence of the two yellow 
curves. 




Figure 2: The graphic for the <f>\ (green), </>2 (blue) and ^3 (red) as function 
of x 3 for a choice of the parameters A\ = 0.2, A2 = 0.8, ipi = 0.9, B = 
1, A = l,.Bi = l,si = — 1,S2 = —1 • Accordingly, the wave vectors for the 
phase and pseudo-phase have been chosen to be a M = (0, 0, 0, 0.33541) and 
= (1.49638, 1, 0, -1.49638), respectively 
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Figure 3: The graphic for the <j>x (green), </>2(blue) and ^3 (red) for a choice 
of the parameters A x = 0.2, A 2 = 0.99, i/n = 20.01, B = 1, A = l,B x = 
1, si = — 1,S2 = — 1. The wave vectors are a M = (0,0,0,-1.58153) and 
^ = (-1.04879, 1,0, 1.04879), respectively 




Figure 4: Projection on the plane (<^i,<fo) of a sample of about 5000 
consecutive values of the field (f> along the axes (0,0,x 3 ), for the same choice 
of parameters as in Figure 2 
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Figure 5: A sample of 200 consecutive spin configurations along the x 3 axis, 
with parameters as in Figure 2. 
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